Cronin Effect in Hadron Production off Nuclei 
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Recent data from RHIC for high-pT hadrons in gold-gold collisions raised again the long standing 
problem of quantitatively understanding the Cronin effect, i.e. nuclear enhancement of high-pT 
hadrons due to multiple interactions in nuclear matter. In nucleus-nucleus collisions this effect 
has to be reliably calculated as baseline for a signal of new physics in high-pT hadron production. 
The only possibility to test models is to compare with available data for pA collisions, however, 
all existing models for the Cronin effect rely on a fit to the data to be explained. We develop a 
phenomenological description based on the light-cone QCD-dipole approach which allows to explain 
available data without fitting to them and to provide predictions for pA collisions at RHIC and LHC. 
We point out that the mechanism causing Cronin effect drastically changes between the energies of 
fixed target experiments and RHIC-LHC. High-pr hadrons are produced incoherently on different 
nucleons at low energies, whereas the production amplitudes interfere if the energy is sufficiently 
high. 
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It was first observed back in 1975 0] that high-pT 
hadrons are not suppressed in proton- nucleus coUisions, 
but produced copiously. This effect named after James 
Cronin demonstrates that bound nucleons cooperate pro- 
ducing high-pT particles. Indeed, it has been soon real- 
ized that multiple interactions which have a steeper than 
linear A-dependence lead to the observed enhancement. 
An adequate interpretation of the Cronin effect has be- 
come especially important recently in connection with 
data from RHIC for high-p^ hadron production in heavy 
ion collisions |^ . The observed suppression factor can 
be understood as a product of two terms. One is due to 
multiple interactions within the colliding nuclei, analo- 
gous to the Cronin effect. The second factor arises from 
final state interaction with the produced medium, the 
properties of which are thus probed. This second factor, 
the main goal of the experiment, can be extracted from 
data only provided that the Cronin effect for nuclear col- 
lisions can be reliably predicted. However, in spite of the 
qualitative understanding of the underlying dynamics of 
this effect, no satisfactory quantitative explanation of ex- 
isting pA data has been suggested so far. Available mod- 
els contain parameters fitted to the data to be explained 
(e.g. see 0, ||, and miss important physics. In this 
paper we suggest a comprehensive description of the dy- 
namics behind the Cronin effect resulting in parameter- 
free predictions which agree with available data. 

First of all, the mechanism of multiple interactions sig- 
nificantly changes with energy. At low energies a high- 
kx parton is produced off different nucleons incoherently, 
while at high energies it becomes a coherent process. This 
is controlled by the coherence length 
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where is the transverse momentum of the parton pro- 
duced at mid rapidity and then hadronizing into the de- 
tected hadron with transverse momentum px- 

For a coherence length which is shorter than the typical 
internuclcon separation, the projectile interacts incoher- 
ently with individual nucleons, just as for e.g. pp scatter- 
ing. However, QCD factorization is violated by multiple 
scattering as discussed e.g. in Q]. Therefore, broadening 
of transverse momentum caused by initial/final interac- 
tions, should not be translated into a modification of the 
parton distribution of the nucleus if the coherence length 
is short. In the opposite limit, i.e. if the coherence length 
is longer than the nuclear radius Ra, factorization ap- 
plies. All amplitudes interfere coherently and result in a 
collective parton distribution of the nucleus. This differ- 
ence is present in all of the various manners in which such 
interactions are discussed. It is e.g. adequate to view a 
nucleus in the nucleus momentum frame as a cloud of 
partons. Those with small x overlap and are no longer 
associated with any individual nucleon. Small x corre- 
sponds to a long Ic ^ l/(a;r7ijv). Again, factorization 
applies, but the nuclear parton distribution is modified. 
The mean transverse momentum of gluons increase |^ 
since their density saturates at small kx |l^ . 

Short coherence length. Broadening of transverse mo- 
mentum of a projectile parton propagating through a nu- 
clear medium is quite a complicated process involving 
rescatterings of the parton accompanied by gluon radia- 
tion. Apparently, this process involves soft interactions 
and cannot be calculated perturbatively. Instead, one 
should rely on phenomenology. Corresponding calcula- 
tions have been performed in in the framework of the 
light-cone QCD dipole approach. The transverse momen- 
tum distribution of partons after propagation through 
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nuclear matter of thickness = J^^dz pa{z) (the 

nuclear density integrated along the parton trajectory at 
impact parameter b) has the form [ pi] , 



dVi d\2 e''^^ ^'''-^^^ n1^{ri,r2) 



X e-i<,(^-i-*='2.x)TA(fc) 



(2) 



Here r2^„(fi, ^2) is the density matrix describing the im- 
pact parameter distribution of the quark in the incident 
hadron, 



(3) 



where (fcg) is the mean value of the parton primordial 
transverse momentum squared. 

The central ingredient of Eq. (||) is the phcnomcnologi- 
cal cross section cr|^(?'T, x) for the interaction of a nucleon 
with a qq dipole of transverse separation tt at Bjorken x. 
In what follows we use the simple parametrization |l2[ 



crqq(rT,x) 



1 



(4) 



the parameters of which were fixed by DIS data: Qg (x) = 
IGeV X (xo/x)^/2 and cto = 23.03 mb; A = 0.288; xq = 
3.04- 10"''. 

Note that the kr distribution of quarks from a single 
q — N scattering process is not singular at kx 0, but 
according to has a Gaussian shape. The phenomenon 



of saturation for soft gluons ||9|, [13| is the driving idea 
of parametrization . Therefore, the mean momentum 
transfer in each scattering is not small, but of the order 
of the saturation scale Qs{x). 

Of course, for projectile gluons the broadening is 
stronger than for quarks and the dipole cross section 
Eq. (0) should be replaced by the glue-glue one CTqq = 

Besides broadening of transverse momentum, initial 
state interactions also lead to energy loss ^A, 15 . While 
induced energy loss in cold nuclear medium is negligibly 
small 
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|] , energy loss due to hadronization in inelas- 
tic scattering reactions (which is basically the same as 
for hadronization in vacuum) is important. The first in- 
elastic interaction of the incident hadron triggers energy 
loss and the parton participating in the high-pT process 
arrives with a noticeably reduced energy |jlj, We 
fixed the energy loss AE to a mean value corresponding 
to the mean path length calculated in ||l^ and a rate of 
energy loss dE/dz — — 2.5GeV/fm. 

For the cross section of pA — > hX at high px we use 
the standard convolution expression based on QCD fac- 
torization iITtIi, 



Ja^RA 
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ij — ^kl 
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h/k 



(5) 



where Fi/p and Fj/A are the distributions of parton 
species i,j in Bjorken xi^2 and transverse momentum in 
the colliding proton and nucleus respectively. However, 
to describe the nonfactorizable multiple interactions the 
beam parton distribution Ff is modified by by the trans- 
verse momentum broadening Eq. (H) and by shifting xx 
to ii = a;i -I- AE/{xiEp). For (^) in (^ we use the 
next-to-leading value from ||^ fitted to data for hadron 
production in pp collisions. For the parton distribution 
functions in a nucleon we use the leading order GRV 
parametrization [ p^ . The nuclear parton distribution, 
Fj/A: is unchanged compared to a free nucleon, except 
at large X2 where it is subject to medium modifications 
(EMC effect) which are parametrized according to [p^ . 
For the hard parton scattering cross section [|l7| we use 
regularization masses ma — 0.8 GeV and niq = 0.2 GcV 
for gluon and quark propagators respectively. Such a 
large effective gluon mass was introduced to reproduce 
the strong nonperturbative light-cone gluon interaction 
|pO| dictated by diffraction data. The fragmentation 
functions of a parton k into the final hadron h, -0^/^ 
are taken from |2^ in leading order. We use the realistic 
Woods-Saxon parametrization for the nuclear density. 

As far as all the parameters in (^) are fitted to data for 
proton target, we have no further adjustable parameters 
and can predict nuclear effect. The results of parameter- 
free calculations for the production of charged pions are 
compared in Fig. Qwith fixed target data. Rw/Be{PT) is 
the ratio of the tungsten and beryllium cross sections at 
200 - 400 GeV pi and 800 GeV l2| as function of pT- 
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FIG. 1: Ratio of the charged pion production cross sections 
for tungsten and beryllium function of the transverse mo- 
mentum of the produced pions. The curves correspond to the 
parameter- free calculation Eq. (^, the data are from fixed 
target experiments |22, M] 



Long coherence length. In the limit of Ic S> Ra a hard 
fluctuation in the incident proton containing a high-p^ 
parton propagates through the whole nucleus and may 
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be freed by the interaction. Since multiple interactions 
in the nucleus supply a larger momentum transfer than 
a nucleon target, they are able to resolve harder fluctua- 
tions, i.e. the average transverse momentum of produced 
hadrons increases. In this case broadening looks like color 
filtering rather than Brownian motion. 

Instead of QCD factorization we employ the light-cone 
dipole formalism in the rest frame of the target which 
leads to another factorized expression, valid at 2:2 ^ 1, 



(pt) 



'(7{GA^ G1G2X)® D, 



h/Gi 
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We assume that high-pT hadrons originate mainly from 
radiated gluons at such high energies. The cross section 
of gluon radiation reads M, Eq, EO], 



da{GA -> G1G2X) 
(FpT dyi 



X *GG(fl,a)^'GG(?^2,a) 



1 „ g-^T,fG(ri.x)TA(fc) 



(7) 



Here a = p+(G'i)/p+(G') is the momentum fraction of 
the radiated gluon; (t|q (r, a) is the dipole cross section 
for a three-gluon colorless system, where r is the trans- 
verse separation of the final gluons Gi and G2 ■ It can be 
expressed in terms of the usual qq dipole cross sections, 



\{^qqi^) + ^qqio^r) + aqq[{l - a)r]| . (8) 



The light-cone wave function of the Gi — G2 Fock com- 
ponent of the incoming gluon including the nonperturba- 
tive interaction of the gluons reads p3 , 



^GG{r,a) 



2rl 



a(er •e)(e2 ■ r) 



+ (1 - a)(e2* • e)(ei ■ r) - a{l - a){el ■ e2*)(e • f) 



(9) 



where rg = 0.3 fm is the parameter characterizing the 
strength of the nonperturbative interaction which was 
fitted to data on diffractive pp scattering. The product 
of the wave functions is averaged in ^ over the initial 
gluon polarization, e, and summed over the final ones, 

6*1, 2- 

Expression with the exponentials expanded to first 
order in the nuclear thickness also provides the cross sec- 
tion for gluon radiation in pp collisions. This cross section 
reproduces well the measured pion spectra in pp colli- 
sions. The results for the ratio of pion production rates 
in pA and pp collisions obtained using Eqs. for 
mid rapidity at the energy of LHC, -y/s = 5.5 TeV are 
shown by curve in Fig. ^ 

Note that at the high LHC energy the eikonal for- 
mula Eq. d?]) is not exact. The higher Fock components 
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FIG. 2: Ratio of p — Au to pp cross sections as function 
of transverse momentum of produced pions at the energy 
of LHC calculated with Eq. (^). The dashed and solid 
curves correspond to calculations without and with gluon 
shadowing respectively. 



jSG), 140), etc. lead to additional corrections called 
gluon shadowing. These fiuctuations are heavier than 
120), correspondingly, the coherence length is shorter, 
and one should sum over all different trajectories of the 
gluons. This problem was solved in HJ, and a 
suppression factor Rg(x,Q^ ,b) due to gluon shadowing 
was derived. Here we make use of those results replacing 
the dipole cross sections in (Q), (T3G by Rgo'3G- This 
suppression factor leads to a reduction of the Cronin ef- 
fect as is demonstrated by the solid curve in Fig. ||. Note 
that this curve approaches unity from below at high px ■ 

Predictions for RHIC. The calculations in the energy 
range of RHIC are most complicated since this is the 
transition region between the regimes of long (small pt) 
and short (large pr) coherence lengths. One can deal 
with this situation relying on the light-cone Green func- 
tion formalism . However, in this case the inte- 
grations involved become too complicated. Fortunately, 
the coherence length at the energy of RHIC is rather 
long, Zc ~ 5 fm, within the pr-range where the Cronin 
effect has an appreciable magnitude. Therefore, the cor- 
rections to the asymptotic expression Eq. (^) should not 
be large and can be approximated by linear interpolation 
performed by means of the the so called nuclear longitu- 
dinal formfactor Fa((7c, 6) H 0, 



aj,A{PT) = j d%[[l-{Fl(q,,b)) 

+ {Fl{q,Malf^-{pT,b)} 



{PT,b) 



(10) 



Here apA(j>T,b) is the unintegrated ^-dependent contri- 
bution to the cross section OpAipr), 



FA{qc,b) 



TA{h) 



dz pA{b, z) e''^'^'' 



(11) 
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where Qc = l/h- The formfactor is averaged weighted 
with the cross section at fixed pt and varying initial and 
final parton momenta. 

Expression (|l^) interpolates between the cross sections 

which are 



'pA (pt), Eq. §) , and a;'=/«-* (pr), Eq. ^ 
shown in Fig. || by dotted and dashed curves respectively. 
It is interesting that the dashed curve exposes a weaker 
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FIG. 3: Predictions for RHIC. The dotted and dashed 
curves are calculated at ^/s = 200 GeV using Eqs. (P) and 
(^) respectively. The final prediction taking into account 
the coherence length is shown by the solid curve. 



nuclear enhancement than the dotted one. This might 
be interpreted as Landau-Pomeranchuk suppression of 
the radiation spectrum compared to the Bethe-Heitler 
regime. 

Our prediction for y/s — 200 GeV calculated with 
Eq. (|l^) is depicted by the solid curve which nearly co- 
incides with the Ic ^ Ra one at pr < 2 GeV and is 
rather close to it at higher px- Ic ^ Ra regime at 
higher pT- Eventually, all three curves approach 1 at 
large px > 10 GeV. 

No sizeable gluon shadowing is expected at RHIC en- 
ergy. The reason is that the effective coherence length 
for gluon shadowing evaluated in ||29|] is nearly an order 
of magnitude shorter than 1^ for single gluon radiation as 
given by (|l]). 

Summary: the mechanism of high-pT hadron produc- 
tion has two limiting regimes. At Ic ^ Ra a high-pT par- 
ticle is produced incoherently on different nucleons, and 
the Cronin effect is due to soft multiple initial/final state 
interactions which break QCD factorization. On the con- 
trary, for Ic ^ Ra the process of gluon radiation takes 
long time even for high transverse momenta. As a result, 
coherent radiation from different nucleons is subject to 
Landau-Pomeranchuk suppression. Using the light-cone 
dipole approach we provided the first parameter-free cal- 
culations for the Cronin effect in pA collisions, i.e. no fit 
is done to the data to be described. Our results agree 
well with available data and we provided predictions for 
high-pT pion production at RHIC and LHC. 
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